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V ■ Abstract 

CO 

We show that the low-energy physics of the spin-tube model in presence of a 
critical magnetic field can be described by a broken SU(3) spin chain. Using 
the Lieb-Schultz-Mattis Theorem we characterize the possible magnetization 
£^ ■ plateaus and study the critical behavior in the region of transition between 

the plateaus (S z ) = 1/2 and (S z ) = 3/2 by means of renormalization group 
calculations performed on the bosonized effective continuum field theory. We 
| show that in certain regions of the parameter space of the effective theory the 

system remains gapless, and we compute the spin-spin correlation functions 
in these regions. We also discuss the possibility of a plateau at (S z ) = 1, and 
show that although there exists in the continuum theory a term that might 
cause the appearance of a plateau there, such term is unlikely to be relevant. 
' This conjecture is proved by DMRG techniques. The modifications of the 

three-leg ladder Hamiltonian that might show plateaus at (S z ) = 1,5/6,7/6 
are discussed, and we give the expected form of correlation functions on the 
g ! (S z ) = 1 plateau. 
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I. INTRODUCTION 



One-dimensional and quasi-one dimensional quantum spin systems have attracted much 
attention in recent years due to the large number of experimental realizations of such systems 
and the variety of theoretical techniques, both analytical and numerical, available to study 
the relevant models. Due to the presence of large quantum fluctuations in low dimensions, 
these systems present unusual properties such as a gap between a singlet ground state 
and excited non-singlet states. Examples include spin ladder systems in which a small 
number of one-dimensional spin-1/2 chains interact among themselvesB. In this case, in a 
way very similar to the Haldane spin-S problem!!, it has been found that if the number of 
chains is even the system effectively behaves as an integer spin chain with a gap in the 
low-energy spectrum, while it remains massless for an odd number of chains. Some two- 
chain ladders which exhibit a gap are SrCv,20$ and Ci^Cs-f^A^) 2C1M, and an example 
of a gapless three-chain ladder is Sr 2 Cu^O^. Thus far we implicitly assumed that the 
boundary conditions in the transverse direction were open boundary conditions (OBC). 
These boundary conditions correspond to having all the chains lying in the same plane. 
This is the situation encountered in experimental systems such as Sr2Cu^O^. In contrast 
with OBC, periodic boundary conditions (PBC) are frustrating for (2n+l) coupled chains. 
As a consequence all the spin excitations are gappedii in the case of periodic boundary 
conditions. They are also gapped for 2n chains with PBC but the mechanism is related to 
singlet formation as in the OBC case and not frustration. The PBC could be achieved in 
an experimental system by having the coupled chains forming a cylinder instead of lying in 
a plane. 

A richer behaviour emerges when these gapped or ungapped systems are placed in a mag- 
netic field. Then it is possible for an integer spin chain to be gapless and a half-odd-integer 
spin chain to show a gap above the ground state for certain values of the fieldQll2l. This 
has been demonstrated by several methods such as bosonizatiqni'0, perturbation theoryE! 
or density-matrix renormalization group method (DMRG)liat3. In particular, it has been 
shown that spin-1/2 chains and ladders with a gap undergo a continuous phase transition 
from a commensurate zero uniform magnetization phase to an incommensurate phase with 
non-zero magnetization^, and the magnetization of the system can exhibit plateaus at cer- 
tain non-zero values of the magnetic fieldiS Further, a striking property of the quantum 
spin-chains in a uniform magnetic field pointing along the direction of the axial symme- 
try (^-direction), is the topological quantization of the magnetization under a changing 
of the magnetic fielJ. It was shown that starting from a generalized Lieb-Schultz-Mattis 
(LSM) theorem0, that translationally invariant spin chains in an applied field can be gapful 
without breaking translational symmetry, only when the magnetization per spin, m, obeys 
S — m = integer, where S is the maximum possible spin in each unit cell of the Hamiltonian. 
Such gapped phases correspond to plateaus at these quantized values of m. In Ref. |TB| the 
behavior of the magnetization versus magnetic field has been investigated in details using 
DMRG techniques for three coupled spin 1/2 chains with both periodic and open boundary 
conditions. Plateaus have been obtained at m = 1/2 and m = 3/2 in agreement with Ref. 
[7|. Further, for the case of PBC a small plateau at m = was also obtainediSH. Finally, 
there seems to exist some weak evidence for a plateau at m = 1 for PBCEI. Strong coupling 
Low Energy Hamiltonians for these two systems were also derived in Ref. [18. 
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In this paper, we investigate a three-chain ladder with PBC (spin-tube) in presence of 
a uniform magnetic field by using bosonization and renormalization group techniques. We 
are concerned by the transition region between the magnetization plateaus at m = 1/2 and 
m = 3/2. Our analysis is based on the low-energy effective Hamiltonian (LEH) derived for 
strong coupling between the rungs@. We identify the LEH as an anisotropic SU(3) spin 
chain with symmetry breaking terms in a longitudinal magnetic field, and analyze its low- 
energy physics via bosonization and RG techniques. This approach allows us to predict 
the behavior of the spin-spin correlation functions in this transition region and the NMR 
relaxation rate. This also allows for an investigation of the possibility of a plateau at m — 1. 

The paper is organized as follows. In Sec.[II], we recall the derivational of the LEH 
and reduce it to an anisotropic SU(3) spin chain, while that in OBC ladder reduces to an 
anisotropic SU(2) spin chain. The difference between PBC and OBC models becomes obvi- 
ous in the language of effective Hamiltonians. In Sec . [HI Cl we discuss the generalization of 
the LSM theoremO for SU(3) spin chains. We then review the analysis of the magnetization 
process of isotropic SU(3) spin chain, and discuss the difficulty^ in obtaining the position of 
the plateaus from the LEH. We show there is no gap at m = 1 in a simple numerical analysis 
using density matrix renormalization group method as well. The bosonized Hamiltonian is 
derived in Sec.|TV|. In section M, we analyze the low-energy effective Hamiltonian in a weak- 
coupling limit by calculating the one loop renormalization group (RG) in the marginally 
perturbed SU(3) Wess-Zumino-Witten mo del@ and discuss the renormalization group flow. 
At weak coupling, the flow is to an invariant surface, leading to gapless excitations above 
the ground state with no breaking of the discrete symmetry. Based on the weak coupling 
renormalization group analysis and the usual continuity between weak coupling and strong 
coupling in one- dimensional systems, we claim that the spin-tube is described by a two 
component Luttinger liquid at low energy and long wavelength. In section[VT|, we discuss 
the effect of a variation of the magnetic field in that problem, and show that it does not 
affect the two component Luttinger liquid behavior. We discuss the analogy of this two 
component Luttinger liquid with the S2 phase of the bilinear-biquadratic spin 1 chainEI. 
Then, having established the equivalence with the two component Luttinger liquid, we cal- 
culate the spin correlation functions in the critical region and the temperature dependence 
of the NMR longitudinal relaxation rate Tf~ . We also present a theoretical description of 
the plateau at m = 1 in the framework of bosonization. Comparing this description with the 
numerical results of Ref. [18] we conclude that the presence of a plateau at m = 1 is unlikely 
in the spin-tube. We verify our results on the absence of plateaus at m = 1 using DMRG, 
and indicate the modifications of the spin-tube Hamiltonian that could lead to a plateau. 
Sec. |Vll| contains the concluding remarks. Technical details can be found in the appendices. 



II. LOW-ENERGY EFFECTIVE HAMILTONIAN OF THE SPIN-TUBE 

The Hamiltonian of the three-chain ladder with periodic boundary conditions (PBC) in 
the presence of an external magnetic field is, 

N 3 _ N 3 ^ ^ 3 _ 

H = J ^ ^ Si )P Si+i tP + J± ^ ^ Si tP Si tP+ i — ^ Si jP , (2-1) 

i=l p=l i=l p=l i=l p=l 
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where p (resp. i) is a chain (resp. site) index, J is the coupling along the chain, J± the 
transverse coupling and the site (z, 4) is identified with the site (i, 1). The three-chain ladder 
with periodic boundary conditions can be viewed as forming a tube with an equilateral 
triangular cross section (see Fig. [I]). We will refer to this system as a spin-tube. 

In the rest of the paper we shall consider the model for J± ^> J, and the aim of this 
section is to recall briefly the derivation of the low energy HamiltonianEl in this limit. To 
begin with, for J = 0, the system consists of independent rungs. The eight states of a given 
rung fall into a spin-3/2 quadruplet and two spin- 1/2 doublets. In the absence of a magnetic 
field, the spin 3/2 states on a given triangle are all degenerate with energy 3/4 J±. These 
states are: 

|3/2;3/2> = |m> 

|3/2;i/2) = i=[| m> + IUT> + UTT>] 
|3/2;-l/2) = i| at) + I Itl) + I TU>] 

|3/2;-3/2} = UU> (2.2) 

Also, in the absence of a magnetic field and on a given rung, the two spin 1/2 doublets, 
corresponding to the left and right chiralities (-/+), are degenerate with energy — 3/4Jj_. 
These states are: 

|t+) = i| ITT) + il TIT) + J 2 I ttl>] 



i+> = -^0 TID+j'I ITI)+j 2 I III)] 
T-) = ^[| ITT) + j 2 I TIT) + j'I TTI)] 
l_> * OTU>+/im>+ilUT}] (2.3) 



where j = exp(-| i ). 

When an external magnetic field is switched on the degeneracy in the different multiplets 
is lifted. The energy levels of the state | |||) (in the spin-3/2 multiplet) and the spin-1/2 
states | t +), | T ~~ ) cross at h c = | J± (see. Fig. |2|). As a result for h < h c one has a ground 
state magnetization (S z ) = 1/2 and for h > h c , (S z ) = 3/2, i.e. h c is a transition point 
between two magnetization plateaus. If a small coupling J is turned on, this transition is 
expected to broaden between /ii/2,+ and ^3/2,-, where h 3 / 2 _ — ^1/2,+ is of the order of J. We 
expect that in this interval (S z ) increases continuously with h. In this limit the properties 
of the system can be studied by perturbing with Hi around the decoupled rung hamiltonian 
Ho, 

H = H + Hi (2.4) 

JV 3 N 3 

Hp — j± Si tP Si )P+ i — h c ^2^2 st, P i (2-5) 

i=l p=l i=l p=l 
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IV 3 N 3 

#i = JEE S i>p S i+1>p - - E E S?, p . (2.6) 

i=l p=l i=l p=l 

At h = h c the groundstate of H is 3^ fold degenerate, the states | j — ) i5 | f 
tTT)i) t° be denoted respectively |2) i; |3);, spanning the low-energy subspace. H 1 
lifts the degeneracy in the subspace, leading to an effective Hamiltonian that can be derived 
by standard perturbation theory. Since in the truncated subspace there are 3 states per 
triangle, it is natural to express the spin operators in the basis given by Gell-Mann matrices 
X a , a = 1 ... 8. (The conventions we use for the matrices can be found for instance in Refs. 



21] and |22|). By considering the action of the spin operators and S* 23 on each state of 
the truncated Hilbert space the spin operators can be expressed in terms of the matrices as, 

S+ P = ^Ij^ + ^+J 2 ^^-*^)} (2-7) 
Sip = i[§J - ^ - /^(A! + iA?) - j^iXl - (2-8) 
where I is the identity matrix. The total rung spin is given by: 

s ' = (I 1 ~ I) • (2 - 9) 

The effective Hamiltonian to first order then becomes: 

H ee . = F + #/, (2.10) 

j N 8 

#o= 7 EE A?A? +1 (2.11) 

4 j=l a =l 

N N h - - N 

I bp 



ieff. 

\ \ + l + " 

i=l i=l i=l 



Hi = ff Am + A 2 A 2 +1 ] + « E Af A* +1 + «'E A^ + i + ^ E A? (2-12) 

V3 i=1 



In our case, g = —J/12, u = —J/4, -u' = —5 J/36 and /i e ff. = h — h c — 5 J/9; hereafter we 
choose our units so that J = 1. The Hamiltonian ( |2.10| ) is written as an isotropic SU(3) 
spin chain H and terms Hj that break the symmetry. This form will be convenient later 
on when we study such questions as what regions of parameter space are gapless and the 
behavior of correlations functions there. 

Another form of the Hamiltonian is convenient when one wishes to study the plateau 
structure. Introduced a different basis of SU(3) operators T-j 4 , T 2 ± , T 3 ± , and T z defined by: 

Tt = (Ax ± zA 2 )/2 (2.13) 
T£ = (A 4 ±*A 5 )/2 (2.14) 
T 3 ± = (A 6 ±2A 7 )/2 (2.15) 

T z = -2-^L (2.16) 
Then, to first order, and up to a constant, the effective Hamiltonian reads: 
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Heff ~ 2 + ^,2^+1,2 + ^3^i+l,3 + ^£,3^+1,3] + 

i 



This is the Hamiltonian derived in Ref. 18. In this form the underlying structure of an 



anisotropic SU(3) spin chain in a "A magnetic field" is unexploited. The correspondance 
between our notations and those of Ref. |H] can be found in table |I[ 

The form of the Hamiltonian (|2.10 ) may help in relating our model to integrable versions 



of the SU(3) spin chains. Isotropic spin chains are known to be integrable by Bethe Ansatz 
techniquescjo. The magnetization process of SU(3) spin chains with a magnetic field cou- 
pled to A 3 has been analyzed in the context of the bilinear-biquadratic spin 1 chain at the 
integrable Uimin-Lai-Sutherland point&H by solving numerically the Bethe Ansatz equa- 
tions. There exist also integrable anisotropic SU(3) spin-chainsca, but the chain described 
by the Hamiltonian ( 2.10Q is not one of them. To study the magnetization process and the 



correlation functions of the chain we will therefore have to resort to a combination of approx- 
imate methods such as bosonization, renormalization group and strong coupling analysis. 
We will, first, analyze the strong-coupling effective Hamiltonian (|2.17 ) to study the region 



of transitions between plateaus of the magnetization, next we will discuss the low-energy 
properties of the effective Hamiltonian ( p.lOj) via a renormalization group analysis. 

We conclude this section by contrasting the Open and Periodic boundary conditions. 
The same strong coupling analysis can be done for the OBC case. In contrast with the PBC 
case, we have only a two fold degeneracy instead of three at J = under a strong field 
h = These two low energy states are 

|3/2;3/2) = |TTT) 

|l/2; 1/2) = J 2 - (\\ TT1) + \\ ITT) - | TIT)-) (2.18) 



The effective Hamiltonian to first order perturbation in J / Jj_ becomes the well known spin- 
1/2 XXZ model, 

Heff = j E + ofrk + Aofo^) -g-^W at, (2.19) 

where A = y|. It is well known that this Hamiltonian has no gap for A < 1 except at the 
saturated magnetization a z = ±1, which corresponds to m — 1/2, 3/2 in the original ladder 
model. This agrees with a weak coupling analysis (J >> Jj_) based on bosonizationill. In 
the strong coupling analysis, one can clearly see the difference between PBC and OBC in 
their effective Hamiltonian. 



III. ANALYSIS OF THE MAGNETIZATION PROCESS OF AN ANISOTROPIC 

SU(3) CHAIN 

We begin the analysis by discussing the generalization of the Lieb, Schultz, Mattis the- 
orem to the SU(3) case@0. The theorem allows us to predict the possible locations of the 
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plateaus. Comparing then these predictions for a chain of SU(3) spins for the plateaus' 
location with those for 3 spin-1/2 chains we conclude that the mapping onto a SU(3) spin 
chain does not introduce spurious plateaus, and should therefore give physically correct re- 
sults for \h — h c \ ~ J. Next, we will discuss the magnetization process of an isotropic SU(3) 
spin chain and address the question whether a cusp appears in the magnetization versus 
magnetic field curve (such a cusp is not related to magnetization plateaus). A cusp has been 
observed in spin-1 chains with bilinear-biquadratic exchange which can be mapped onto an 
SU(3) spin chain for a special value of the biquadratic exchange. We argue that in our case 
a cusp in the magnetization should not be expected. Finally, We recall how the values of 
the magnetic field corresponding to the plateaus^l at (S z ) = 1/2, 3/2 are obtained from the 
anisotropic SU(3) spin chain. 

A. The Lieb-Schultz-Mattis Theorem for a chain of SU(3) spins 

Consider an anisotropic chain of SU(3) spins, 



N 



(3.1) 



i=l L a 

where CL\ — Qj2i CI4. — O5, Clg — O7, which we rewrite as, 

H = 2£ [ai(7?i^+i,i + + 04(^3^1,2 + 1^^1,2) + a6(^+i,3 + ^+1,3)+ 

i 

+ 2a 3 AfA| fl + 2a 8 A^ +1 + 2/?(A, 3 A* +1 + AfAf +1 )] . (3.2) 

The anisotropic SU(3) chain we are considering in this paper falls into this class of Hamil- 
tonians. 

The purpose of the LSM theorem is to classify possible gapless excitations above the 
ground state. Introduce the operators: 

/ cs N \ 

U 3 = exp \ i— 5^ nAfj (3.3) 

^= exp ( ! ^sS" A ") (3 ' 4) 

and begin by studying the state U 3 \0). We wish to compare its energy expectation value 
with the vacuum's. Consider therefore the expression (0\U j 3 HU 3 - H\0). Using, 

+ 1 o„ Inn 1 

UlT+ 2 U 3 = e^T+ 2 

UlT+ 3 U 3 = (3.5) 

we find upon developing it as a power series in 1/A (N is the size of the system), that the 
zeroth order term vanishes, since it contains averages of the form, i{T^ +ltP T~ — T~ +1)P T+) , 
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which cancel when invariance under parity m — > 2n + 1 — m, is used: (Tn+i, p T np ) = 
(^n+i^np)' Therefore, the expansion begins with a first order term in 1/N, 

^ x (~^T £ [ 4 «i (^+1,1^1 + H.c.) + a 2 (T+ +1>2 T- 2 + if.c.) + a^T^T-g + F.c.)] j 

Translation invariance guarantees that the coefficient of 1/N in this term is indeed finite. 
The translation operator T, is defined by: 

T- 1 KT={ X ^ (36) 



and translation invariance, T|0) = |0), (together with the expression of C/3, Eq.( |3.3[ )), implies: 
TC/ 3 |0) = TUzT^TlO) = e - 2 "( A ^- m 3)|o), where m 3 = ^En=i A n- Therefore, the state 
Us I 0) is orthogonal to the ground state if m 3 is non-integer. This state has energy 0(1/N) 
above the ground state, indicating either a ground state with a broken symmetry or gapless 
excitations above the ground state for m non-intege 101. A gap in the excitation spectrum 
can only exist for m 3 = —1,0, 1 in the absence of broken symmetry ground states. For a 
ground state with p-site-periodicity instead of one site translational symmetry, 

T*|0) = |0>, (3.7) 

we obtain, T p [/ 3 |0) = e~ 2m( ^^o ^-i - *" 713 ) f or a p OS jtive integer p. A gap excitation can then 
appear only for m 3 = q/p with an integer q = —p, ■ ■ ■ ,p. 
If we now consider the action of Us, we have: 

U 8 T^ }1 U 8 = T+i 
UlT+ 2 U 8 = e^T+ 2 



UlT+ 3 U 8 = e^T+ 3 (3.8) 

Once again, (0 | U ] 8 HU 8 - H | 0) = 0(1/N), but this time: TU 8 | 0) = e~ l ^ {X *~ m&) U 8 | 0), 
where m 8 = Yln=i ^n- This implies that a gapful excitation on a translational invariant 
ground state is only possible for ^| = n + 1/3, where n is integer. As in the case of U 3 , 
a gap on p-site-periodic ground state can exist for = n + 1/3. As a consequence, 
two conditions have to be met to avoid having gapless excitations above the translationally 
invariant ground state (p=l): 

7718 , 1 /q 

7i = " +1/3 

m 3 = n'. 

The preceding discussion is quite general. In the specific case we are considering, the mag- 
netic excitations are described by S z = 5/6 — A 8 /\/3. Therefore, we expect that the magne- 
tization plateaus are associated with the absence of excited states above the ground states 
generated by U 8 . This implies that the only possible magnetization plateaus correspond to 
A 8 . Clearly, such a conclusion could have been reached by considering the LSM Theorem 



. 2 7T n 
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for a three chain system. However, the LSM theorem for the SU(3) chain also indicates 
the possibility of a gap for chiral excitations for (A3) = —1,0, 1. Because the Hamiltonian 
( f2.11| )-( |2.12j ) shows chiral symmetry, one has necessarily (A3) = 0. Therefore, according to 
the LSM theorem, we are in a situation where a gap in the chirality modes can be present 
and we should decide whether this gap actually opens. Let us consider two simple limiting 
cases. For (X 8 ) / y/3 = 1/3, i.e. (S z ) = 1/2, the only possible states on each site of the SU(3) 
spin chain are | 1) and| 2). In that case, it can be easily seen that the Hamiltonian fl2.17p 
reduces to an effective XY chain@, implying gapless chiral excitations. On the other hand, 
if (Ag)/\/3 = —2/3, the only remaining mode on each site is |3). One has indeed in this 
state A3 = 0, but this time there is a gap in chiral excitations. If there is a ground state with 
broken translational symmetry, the situation is more complicated. Therefore, the question 
of actual gap opening cannot be settled by the LSM theorem alone. 



B. Comparison with the magnetization process of an isotropic SU(3) spin chain 

The isotropic SU(3) spin chain is known to be integrable by the Bethe Ansatzffl. It is 
also known that the bilinear biquadratic spin-1 chain defined by the Hamiltonian: 

H = £ [Si.Si+i + /3(^4+i) 2 - hSf] (3.9) 

i 

for (3 — 1 (ULS point) can be mapped onto an isotropic SU(3) spin chain. In the context of 
the bilinear-biquadratic spin chain, the magnetization process of isotropic SU(3) spin chains 
has been investigated in detail by solving numerically the Bethe Ansatz equations. In that 
case, the magnetic field couples to A3 = S z . A cusp in was obtained in the magnetization as 
a function of the magnetic field@@. The cusp is due to the fact that three different types 
of Bethe Ansatz quasiparticles have respectively chemical potentials: h, 0, —h. As a result, 
when h is large enough the band with the highest chemical potential is emptied causing the 
cusp in the magnetization. Also the effect of an anisotropy —D(S*) 2 oc A 8 on the bilinear 
biquadratic spin-1 chain at the ULS point was studiedc3. Rephrasing the results in the 
context of SU(3) spin chains, it was found that when one applies a field that couples to A 8 
(respectively A 3 + ^=, A 3 — ^) the resulting average value of A 8 (respectively A 3 + ^?=, A 3 — ^|) 
shows no cusp as a function of the applied magnetic field. The reason is that in this case, the 
Bethe ansatz particles have chemical potentials — h, —h, 2h which prevent band emptying 
effects. In our case, similarly, the magnetic field couples to A 8 . We thus conjecture that 
although the anisotropy renders the system non-integrable, the absence of band emptying 
effects should persist in the anisotropic case preventing any cusp in the magnetization. 
Also, on general grounds, the Hamiltonian ( p,17|) is invariant under exchange of chirali- 



ties. Therefore, we expect to have (A 3 ) = 0, whatever the applied magnetic field. 



C. Strong-coupling analysis of the effective Hamiltonian: Magnetization plateaus 

In this section we will recall the evolution of the ground state magnetization as a function 
of the external magnetic field. As noted before, when the intrachain coupling J is set to 
zero, upon switching the external magnetic field on, we find at increasing h that the ground 
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state of a given rung undergoes a transition between the spin-3/2 state, | fit) an d the 
spin-1/2 states, | | +), | j — ) at /i c = resulting in a sudden jump of the magnetization 
between m = 1/2 and m = 3/2 (as shown in Fig. |2|). If the coupling J is non-zero but small 
this transition is broadened between h 3 / 2 - and ^1/2,+ which can be identified, respectively, 
with the lower and upper critical fields of the saturation plateaus of the magnetization at 



m = 3/2 and m = 1/2 in the terminology of ref. |18[ As noted there, it is easy to obtain 

h /2 ,- = (p± + 2J). (3.10) 

by considering the condition for the ferromagnetic state | 33 • • • 3) to be stable upon the 
introduction of spin 1/2 states |1) or |2). 

It is harder to determine the lower critical field hi/2,+, below which the magnetization 
plateau is at m = 1/2. Indeed, on the plateau m = 1/2, there are two possible states 
on each site, |1) or |2), so that this plateau is described by an effective XY chain. As a 
consequence, the ground state wavefunction for m = 1/2, obtained from the Jordan- Wigner 
transformation is a complicated linear combinations of states of the form | • ■ ■ 22121- ■ ■). 
One has to consider the energy loss created by the introduction of a |3) state in that chain, 
and balance it with the energy gained from the magnetic field. This problem bears some 
similarity to the dynamics of a few holes in a t-J chain, which has an SU(2) xU(l) symmetry 
instead of the U(l) x U(l) symmetry in our model. The analogy with the t-J model suggests 
a two-component Luttinger liquid behavior of the system in a large part of the phase diagram. 



IV. BOSONIZATION AND WEAK-COUPLING ANALYSIS 

We proceed now to study the long distance properties of the effective Hamiltonian H c r 
defined in Eq. ( |2.10| ). It is a sum of an isotropic SU(3) spin chain Hamiltonian plus SU(3) 
symmetry breaking terms, 



with 



H cS . = H + H 1 + H 2 + H 3 + H h (4.1) 



j * 

^o = tE E W +1 (4.2) 

^ j=i o=l,...,8 
N 

ffi = ?EW4i + %l ( 4 - 3 ) 

i=i 

N 

H 2 = uY,^ l+ i (4-4) 

i=l 
N 

ffs^EWi ( 4 - 5 ) 

1=1 

N A 8 

Hh = K ff J2 ^ (4-6) 

In our case, q = —1/12 J, u = — 1/4J, u' = — 5/36J, and h e ff = — -^J; hereafter we 
choose our units so that J — 1. 
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A. Non-abelian bosonization of a SU(3) spin chain 



The SU(3) invariant Hamiltonian Hq can be solved exactly by the Bethe ansat2 
The solution shows that the SU(3) spin chain has 2 branches of excitations, with dispersion 
e i( fc ) = 4 sin(ti/3) [cos(7rj/3 - \k\) - cos(7rj/3)], j = 1,2. These excitations are gapless, 
and for \k\ — > 0, one has ei(k) = 62(h) ~ ^-j \ k |, i.e. the dispersion relation assumes at 
long wavelength a massless relativistic form. Accordingly, the low energy, long wavelength 
excitations of the SU(3) spin chain can be bosonized. More precisely, these excitations 
are described^ by the SU(3) level 1 (£1/(3)1) Wess-Zumino-Novikov-Witten (WZNW) 
modelil, perturbed by a marginally irrelevant SU(3) invariant operator. A review of WZNW 
models can be found in Ref. In Hamiltonian form, the 57/(3) 1 model can be written as: 



2tt 31 

Hwznw = T E : JrWJrW ■ + ■ J a L(x)Jl(x) ■■ (4-7) 

6 a=l 

where the right and left currents satisfy the following commutation relations (Kac-Moody 
algebra at level 1): 

[Jr(l)i J r {L) ] = H^Hx - v)Jl {L) {y) + ^'(x - y). (4.8) 

In Equation (|4.8|) , the f al31 are the structure constants of SU(3). The central charge is 
C = lx fl; 1] = 2, indicating that the SU(3) 1 WZNW model can be described in terms of 
two free bosonic fields. As mentioned above, the 57/ (3) spin chain is described asymptotically 
by the SU(3)\ model perturbed by a marginally irrelevant SU(3) invariant term, 

/dx 
7T®°( X )> ( 4 - 9 ) 
271 

where the marginal operator, $o(x) = J2a=i /i?( :r )/L( a; )? couples the right and left currents. 

The finite size correction to the ground state energy of the SU(3) chain can be obtained 
from the Bethe Ansatz solution. These corrections are logarithmic and are in agreement 
with those obtained from the continuum Hamiltonian ( f4.9| ). This situation is very similar 
to the more familiar case of the SU(2) spin chain, which is described at low energy and 
long wavelength by the marginally perturbed SU(2)% WZNW modefl. In general, the 
magnitude of go cannot be obtained from the lattice Hamiltonian in the case of a SU(N) 
spin chain (see the discussion of the case N = 2 in Ref. ^). This is even more problematic 
when one adds perturbations to the SU(3) invariant spin chain. Another difficulty is that 
these perturbations are not small in our case and strictly speaking cannot be treated in 
perturbation theory. However, in one dimension weak and strong coupling behavior are 
often continuously connected&0 so that a weak coupling analysis can provide very valuable 
information on the qualitative physics at strong coupling. Therefore, if we can find a weak 
coupling model that is described by marginally perturbed SU(3)i WZNW model and if we 
add to it small perturbations of the form ( |4.3| )-( f4/5| ), we will be able to make a reasonable 
guess of the low energy long wavelength continuum theory associated with the Hamiltonian 
i/ e ff .• By analogy with the Heisenberg model, we expect that the difference between the weak 
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and the strong coupling regime will reduce to a renormalization of some parameters of the 
effective low energy theory. For non-integrable models, these parameters can be obtained 
numerically by calculating thermodynamic quantities via exact diagonalization methodsc3~@. 

In our case, it is not difficult to see that the spin sector of the SU (3) Hubbard modelUli! 
is a good candidate for a weak coupling model. This model is defined by the Hamiltonian: 

H = -t (4+l,n C i,n + h.c.) + U n i,nni, m (4.10) 

i,n=l,2,3 i,nj^m 

where Cj )n annihilates a fermion of flavor n £ [1,2,3] on site i, and nj jT1 = c| n Cj )n . The 
basic idea is that, starting from the lattice Hamiltonian of the «S77(3) Hubbard model it is 
possible to take the continuum limit and then separate the spin excitations from the charge 
excitations by means of weak coupling bosonization. In the strong coupling limit, U — > <x>, 
a constraint of one fermion per site is imposed, 

E C U^,n=l. (4.11) 

n 

With one fermion per site, the charge degrees of freedom are frozen out and one is left only 
with SU(3) spin degrees of freedom. 

Second order perturbation theory in t then shows that the model can be mapped onto 
an isotropic SU(3) spin chain with the lattice SU(3) spin operators 

A-i = C i,m^n,m C i,m (4-12) 



under the constraint ( I4.11|) . Under the hypothesis of continuity, the same SU(3)\ field 



theory should describe the weak and strong coupling limits in the spin sector. The difference 
between the weak and strong coupling limit corresponds to the disappearance of the charge 
sector. This reduction of the number of degrees of freedom can be obtained in a consistent 
way by treating the constraint (4.11) within the effective theorjSiH. 



Thus, our strong coupling theory is the spin sector of the SU (3) Hubbard model with a 
filling of one fermion per site and U — > oo. Let us discuss the weak coupling regime. The 
constraint ( I4.11|) sets the Fermi momentum at &f = | for the three fermion flavors. Since 



we are interested in low-energy, long wavelength properties, we linearize the spectrum for 
each flavor around the two Fermi points and introduce the right and left moving fermion 
modes in the continuum limit, 

c[ n ~ Me lkFX ^L, n (x) + e-^4»), (4.13) 

where x = ia and a is the lattice spacing. 
For U = 0, the linearized Hamiltonian is: 

-^linearized = ~IVf ^ / ^(V^rAV^n ~ ^L,n d xtpL,n) ■ (4.14) 

This Hamiltonian is conformal invariant and can be rewritten in terms of the right and left 
charge currents Jr(l) = J2n yml) n ^R{L),ni an d the eight SU(3) spin currents (right and left) 

Jr(l) = Sn%(L)n J 2 !!1, /'R(i),m■ ® ne ^ nus separates the charge and spin sectors: 
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H = -^charge + -^spin (4-15) 

where the charge sector, 

-^charge = v F J dx : J R (x)J R (x) : + : J L (x)J L (x) : (4.16) 
and the spin sector is again described by the S77(3)i model discussed earlier, 

H spin =v F J2 I dx: J a R {x)J a R {x) : + : J a L {x)J a L {x) : (4.17) 

The charge currents satisfy U(l) Kac-Moody algebra, whereas the spin currents satisfy 
the SU Kac-Moody algebra as can be checked explicitly. When the interaction is weakly 
turned on, U/t <C 1, it does not break spin charge separation but induces a go oc U termed. 

We discussed thus far non-abelian bosonization in order to stay close to the literature on 
SU (3) spin chains. However, an abelian bosonization approach to the isotropic SU (3) spin 
chains starting from the SU (3) Hubbard model is perfectly feasible. Such an approach has 
been introduced for isotropic SU(N) spin chains in Ref. 43. It is outlined in appendix 
In fact, for the rest of this section, we shall employ abelian bosonization because it renders 
the calculation of correlation functions extremely easy even when the SU (3) symmetry is 
explicitly broken. 



B. Abelian bosonization approach 



Abelian bosonization gives the following Hamiltonian for an ST/ (3) invariant spin chain 
(or the spin sector of the SU(3) Hubbard model): 



H. 



su(3) 



dx 

'1 

2Ua 



-u 



^n a ) 2 + (7TL) 2 + {d x <\>af + {d x c 



dx(cos V8<f) a + cos V2(<f) a + V3~4>b) + cos V2(<f) a - \/3> 6 )) 



(d x (j)a) 2 + (0 x (j) b 



A derivation can be found in Appendix The free term corresponds to Eq.( 4.17|) . 
Under renormalization, H su ^ flows to a fixed point HamiltonianS 



H I —u* 

2vr 



^nj 2 + (vrn fe ) 2 + (dM 2 + {d x (j> b f 



(4.18) 



(4.19) 



where u* is given by the Bethe Ansatz as u* = ^4. One can check using the expressions 



3 4' 



x 



These 



( |A12| ) that this leads to a scaling dimension of 1 for the uniform component of A a ( 
a _1 A", x = ia (a = 1---8), and 2/3 for the 27r/3 component (see Eq. (|A12 H. 
scaling dimensions coincide with those obtained from non-abelian bosonization0 

Turning now to the SU (3) symmetry breaking terms, we find that in the abelian bosoniza- 
tion representation they take the following form: 
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#1 

H 3 



dx 
dx 
dx 

7i 2 



qa 



qa 



qa 



5ua, n , , 9 ua, n , , 9 1 2ua 
— (d x <P a ) 2 + —(d x <f )b ) 2 + Jr ^ 
2 6 J 4(7ra) / 

u'a,~ , .0 5u'a 2 1 2u'a 
-^(d x (j) a ) 2 + —-(d x (j) b ) 2 
o 2 



+ 



2qa 
(2ixa) 



Au'a 



3(7ra) 



3(27ra) 2 
<ix cos v2(j)a cos + 



dx cos a) + 
dx cos(\/8(J) a ) — 



J dx COS V8(f) a + 

\/~3u 



\^3qa 
(2^)2 



(27T) 



dxd x (f) bj 
(4.20) 



u'y/2 



3tT 2 



Or 



H h = -(h - h c ) / dx — (d x (, 

7T 



V2, 



(4.21) 
(4.22) 



The physical interpretation of the terms proportional to d x (j) b is very simple. The bosonized 
Hamiltonian is derived under the assumption that the magnetization per triangle is close to 
5/(6o). When the magnetization per triangle is exactly 5/ (6a) the terms d x (f) b do not appear 
in the Hamiltonian. Therefore, the presence in the Hamiltonian of such terms means that 
the magnetic field needed to impose a magnetization of 5/ (6a) per triangle is renormalized 
away from its bare value. Also, since the Hamiltonian preserves the symmetry between 
+ and — chiralities, it is invariant under the transformation 7r a — > — 7r a , <p a — > — a . In 
particular, this precludes the terms d x cj) a in the Hamiltonian. 

Assembling all terms we finally have the following field theory describing the spin sector 
of the SU(3) Hubbard model in the presence of symmetry breaking perturbations, 



H = v F J2 J ^[(W) 2 + (<9^) 2 ] + +t|^ J dxcosVsM* 



i=a,b 



+ 



4#2 



(2?ra) 



J dx cos y/2(J) a (x) cosVE(j)b(x) + ^| J dx{d x <f> a ) 2 + ^| J dx{d x cj) b f + — J dxd x (p b (4.23) 



with vp = 2tasin(k F a) = \/3ta, h = (^pw + ^u' + ^g), and t the hopping amplitude. In 
our units, t = 1. 

The notations can be made more compact by introducing the vectors <p — (4>a, <&») an d 
ai = (1, 0) , a 2 = (1/2, y/3/2) and a 3 = (1/2, -V3/2), where, 



K„ 



u„ 



K h 



qa 
2nvp 



1 - 



Ua 



qa 



j u'a j 5ua 



ga 



[/a 



1/2 



ga u'a hua 
+ 7: + 



TTVp KVp 3tIVf 7TVj? 
1/2 



1/2 



C/o j 5n'a j ua _|_ qa 



U b =V F [l 



ttvf ttvf 3nvp 6nvp 

Ua 5u'a ua 



qa 



1/2 



TTVp TTVp 3lTV F 6TTVF 

The Hamiltonian can then be rewritten^ as, 



(4.24) 
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The interactions {d x <p a ^) 2 can render the marginal operators cos(v / 80a), cos(\/20 a ± V64>b) 
marginally relevant and cause the opening of a gap. In such a case, the low energy properties 
of the system cannot be described by two massless bosons. One can have either a massive 
and a massless boson or two massive bosons. This depends on the coupling constants u, u', 
m, q, and the magnetic field h. In order to explore this possibility in more detail, one has 
to use renormalization group equations. This is the subject of the forthcoming sections. 



V. THE RENORMALIZATION GROUP FLOW IN ZERO MAGNETIC FIELD 



In this section, we discuss the flow of the renormalization group equation and the phase 
diagram that results from it. Qualitatively, the renormalizationgroup equations are similar 
to the Kosterlitz-Thouless renormalization group equationsE3@ We expect therefore to 
obtain a gapless phase corresponding to the flow to a fixed hypersurface of the 6-dimensional 
space of coupling constants and one (or possibly many) gapped phase where the coupling 
constants flow to infinity. We also expect that the phase transition is of infinite ordeio. 
Our task is therefore to determine the initial conditions and follow the flow. This will allow 
us to conclude on the nature of the ground state of the anisotropic SU(3) chain. 

A straightforward application to the Hamiltonian ( 4.23| ) of the standard method^ 
would be inconvenient since one needs to expand to third order of correlation functions 
in order to get the full one loop RG equations^. We will use instead operator product 



expansion (OPE) techniques&JIH In our case, the algebra of operators (d x (p a ) 2 , (d x 4>b) 2 and 
close under OPE (for details see Appendix O). In particular, 



COS^VoCtj 

cos(VSa^)(x, t)) cos(\/8c50(O, 0)) 
cos (V8acf(x, r)) cos(v^/30(O, 0)) 



-2a 4 



(X 2 + (UT 




(a p l ) 2 (x 2 (d x <f )p ) 2 + r 2 («9 r p ) 2 + 2xrd x <f )p d T <P P 

p=a,b 

-Ka.0 

cos v / 8(a + /?).0(O,O) + 



cos V8(a-/?).0(O,O) 



lead to the following RG equations (see appendix 

tym — 2/1/2 



dyi 
dl 



^ = (1A/4 + 3/22/5)2/2-2/12/2/2 

dl 



dyi 
dl 



yl/i + vl/* 



-df = 3/%2 ' 



(5.1) 
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where we have denoted yi = gi/irvF, with vf the Fermi velocity. Note that is a fixed surface 
{gii 92) — (0,0), because of three truly marginal operators d x (p a d x 4>b, (d x (f> a ) 2 and (d x 4>b) 2 . 

An alternative approach based on non-abelian bosonization can be used. In this ap- 
proach, having expressed the Hamiltonian in terms of products of right and left moving 
currents an operator product expansion for currents is derived^. Such an approach 

leads to the same RG equations as the Abelian bosonization approach. 

The initial values of the running coupling constants (at the cut off scale a) for the spin 
sector of the SU(3) Hubbard model perturbed by -^1,2,3 are given by, 

( \ #i( Q ) ijj , u ' a , w \ 
?/i (a) = = {Ua + qa — + ua)/iTV F ) 

TTVf o 

1 \ 92(a) u'a 

y 2 (a) = = {Ua + —)ixv F 

tivf 3 

/ \ 9&{°) , Ua 5ua qa u'a. . 

y*( a ) = -1— = (—5- + -7T + v + it ) ™f 
irvp z z z b 

(75(a) Ua ua 5u'a qa 

ysw = — — = (—5- + it + + to)/ 1TVf - ( 5 - 2 ) 

nvp Z Z 12 



In the expression for the initial coupling constants (|5.2|) we have u = —J/4, q = —J/12, 
u' = — 5J/36, and we assume J,U <C t. Hereafter, we choose J = 4 and put Vp equal to 
unity, thus the numerical starting values are given by 

y x = -0.365467 + y 
y 2 = -0.0589463 + y Q 
y 4 = -0.878299 - y /2 

y 5 = -0.5039907 - 2/0/2. (5.3) 

Here j/o ^ U. These values are not small, so that the one loop RG equations are not valid. 
However, numerically solving these RG equations with initial conditions ( |5.3| ) shows that 
they flow to a fixed point on the surface (gi, g 2 ) = (0, 0) for any y G [0, 1] (see Figs. § and 
|J). At this fixed point, one has a renormalized Hamiltonian with: 



H 



dx 
2^ 



77 II 



(5.4) 



This proves that certainly at weak coupling the long distance properties of the system are 
described by a two component Luttinger liquid. At strong coupling, i. e. in the case of the 
spin tube, this is only an indication that the two component Luttinger liquid is possible. 
In order to give a definitive proof, one should prove that there is no singularity in the 
ground state energy as couplings increase. Comparing Fig. |3] and f|, one can see that the 
magnitude of the fixed point values of g± and g§ depends on the strength of the marginal 
SU(3) symmetric interaction y . This fact, combined with the fact that the RG equations are 
only valid at weak coupling precludes the use of the RG to give an accurate estimate of K* 
and K^. However, one can still determine from the RG equations whether these quantities 
are larger or smaller than one. Although we stress that these figures should not be taken too 
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seriously, we find using the expressions of K a and Kb as a function of and g$ K* = 1.9 and 
Kl = 1.5 i.e. both are larger than 1. Concerning the question whether the two component 
Luttinger liquid persists at large coupling, we can remark that the deviation from isotropy 
in our case makes the interaction between the SU(3) spins less antiferromagnetic. It is well 
known that in the case of the XXZ chain, reducing the antiferromagnetic character of the 
spin spin interaction (i.e. working at J z < J ) prevents the formation of a gapEI. Therefore, 
it seems likely that no gap would develop in the spectrum. To test our conjecture, numerical 
work, especially calculation of K* a ,Kl by exact diagonalization would prove very valuable. 

The existence of a two component Luttinger liquid phase has important consequences. 
In particular, it implies a non-zero magnetic susceptibility x ^ K b /ub, and a T linear specific 
heat of the form: 

C = f + f. (5.5) 
bu a bu b 

The calculation of the correlation functions and NMR relaxation rate are postponed to 



section VI 



VI. STRONG MAGNETIC FIELD CASE: FIXED POINT HAMILTONIAN AND 

CORRELATION FUNCTIONS 

A. Generic magnetic field 

1. Renormalization group flow under a magnetic field 

Till now, we haven't taken into account the terms associated with the magnetic field 
hb, which can be treated as a perturbation having fixed the external magnetic field at h c = 
3/2J±. To see if the flow remains unchanged in this case, let us reobtain the renormalization 
group equations with finite h. The simplest way to address this problem is to perform a 
Legendre transformationil on the Hamiltonian ( [4.231 ). The non-zero average value of the 



field <pb due to the finite magnetization can be eliminated by a simple shift of the <f>b fields, 
i.e. <pb = <Pb — TTiribX, where m& = — (d x (pb) / n ■ One has the relation between nib and the 
magnetization : 




The cosine terms, however, are not invariant under this shift and the renormalization 
group equations Q5.1]) for the couplings g\, g-i and #3, for a change of the length scale a — > e l a, 
now become (the details of the calculation can be found in the appendix 0): 

= 23/12/4 - ylJ (irm b (l)aV3)/2 

at 
dy 

—jj- = (l/2t/ 4 + 3/23/5)3/2 - y 'W2 -J r o(Trm b (l) a— ■) /2 
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^ = yl/2 + ylJ (nm b (l)a^)/A 

^ = 3/AylJ (nm b (l)a^-), (6.2) 

where Jo is the Bessel function that results from the use of a sharp cutoff in the real space. 
One also has m&(Z) = m&(0)e . One can check that setting m&(0) = one recovers the 



equations (5J). If the RG equation for the magnetization is trivial, the magnetic field on 



the other hand, satisfies a non-trivial RG equation: 



-jr = h + ^=ylJi{r(VQm b (l)a). (6.3) 
dl av c 



Let us discuss qualitatively the physics predicted by the Eqs. (|6.2j). One sees rather easily 
that for m b (l)a <C 1, the Bessel functions tend to zero, so that one is left with a Sine-Gordon 
renormalization group equation for 2/1,2/4. Compared to the case of zero magnetization, we 
see that y 4 is more negative and y\ is smaller in absolute value. Therefore, we expect that 
yi will be even more irrelevant in the presence of the finite magnetization. We conclude 
that the presence of a non-zero magnetization does not affect the two component Luttinger 
liquid behavior. The crossover scale can be roughly estimated as: 

I* ~ ). (6.4) 

m b a 

At this crossover scale, the flow of y$ is completely cut. This implies a variation of K a , K b 
with the magnetization. 

At the value of I given by ( |6.4|) the magnetic energy is of the order of energy cut-off, 
therefore the magnetic field term cannot be treated as a perturbation. When the initial 
magnetization goes to infinity the renormalization is stopped for smaller and smaller /. The 
coupling constants gi,g2,g3 then become zero, while #4,(75 assume the values they have 
at the scale /*. Returning to the Hamiltonian ( |4.23| ), we see that it becomes a quadratic 
Hamiltonian. 



2. Fixed point Hamiltonian 

Following the the preceding discussion, we conclude that the asymptotic behavior of the 
three chain system under a magnetic field is governed by the Hamiltonian: 

H* = I pv F UnU a ) 2 + {d x <Pa? + (vrn,) 2 + (dj b ) 2 ] + 9 -\ f dx{d x <P a f + 9 -\f dx{d x <p b f 

where g\ 5 are functions of the magnetic field. The field <j) b is related to (f) b in the following 
way: 

5 [3 

b = b + vr(m- — )W-x, (6.5) 
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while the dual fields 6 a and 6 b are not shifted. This condition guarantees that <p b satisfies 
periodic boundary conditions. 

The fixed point Hamiltonian can be rewritten: 




where u*K* = vp] i = a,b and u* a b /K* b = vp+2gl 5 /n. Both the velocities of excitations, u i: 
and the compactification radii, Ki, depend on the magnetic field h through g*{h). Therefore, 
the low energy properties of the system are described by two decoupled c = 1 conformal 
field theories with velocities and compactification radii depending on the applied magnetic 
field. 

This is valid at the level of perturbation theory for the spin sector of the SU(3) Hub- 
bard model. However, we are actually interested in the SU(3) anisotropic spin chain for 
which perturbation theory does not apply. In the latter case, we expect, relying on the 
continuity between the weak and the strong coupling regime, that the anisotropic SU(3) 
spin chain under magnetic field will also be described by a two component Luttinger liquid. 
However, the velocities and compactification radii cannot be obtained by perturbation the- 
ory techniques. Nevertheless, it is known that the velocities and compactification radii can 
be obtained by calculating only thermodynamic quantities using, for instance, exact diag- 
onalization techniquesc3'El. The problem of the determination of these exponents in terms 
of measurable thermodynamic quantities in the specific case of the anisotropic SU(3) spin 
chain will be discussed in Appendix [D]. The knowledge of the exponents then permits the 
calculation of correlation functions. This is the subject of the forthcoming section. 



3. Correlation functions 
In this section, we want to calculate the three Matsubara correlation functions: 

Xzz (x,r) = (T t S'(x,t)S*(0,0)) (6.7) 
X+ _, p (x,r) = (T t S+{x,t)S;(0,Q)) (6.8) 
XzzA^r) = (T T S z p (x,r)S z p (0,0)), (6.9) 

where p — 1, 2, 3 is a chain index. The first correlation function is useful for neutron 
scattering experiments, whereas the correlation functions (|6.8j ) are useful for the calculation 
of NMR relaxation rates. The Matsubara correlation functions in Fourier space are given 
by: 

XiM^T) = f drdxex P ^^(T T [S l (x,r),S^O,0))) T , (6.10) 

•/ 

from which the finite temperature correlations are obtained by the analytic continuation 
iuj n — > ui + 20+ . We will first concentrate on the T = calculation, then explain how to 
extend the calculation to finite temperature. 

We begin with the calculation of Xzz- Using the equation ( |2.9| ) we have: 
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Xzz = i(T T A 8 (x ) r)A 8 (0,0)) + ((^)) 2 . 



(6.11) 



Using the bosonized expressions of the SU(3) spins, Eq. ( |A12j ), and the usual expression 
for the bosonized correlation functions^, we obtainEl: 



Ku (v^t) 2 - r 2 P<^-< m --k)) x 



Kg 

6 



37T 2 (x 2 + (u b T) 2 ) 2 



+ 



3(7ra) 



x 2 + (u b ry 



6(7ra) 2 \x 2 + (u a r) 2 J \x 2 + (u b r) 2 
+ c.c, 



2Kj 
3 



(6.12) 



where m = (S z )/a. 

Turning to X+-,pi it is easily seen using Eq. ( |2.8| ) that it is independent of p and equal 

to: 



X+-, p (x,t) = — 



(T T (A X + tA 2 )(x, r)(A x - *A 2 )(0, 0)) + (T T (A 4 + tA 5 )(x, r)(A 4 - *A 5 )(0, 0)) 



Similarly, Xzz, P is independent of p (see Eq. (|2.8| ) and the contribution not already included 
in Xzz is of the form: 



(T T (A 6 + iA 7 )(x, r)(A 6 - iA 7 )(0, 0)). 
The expressions of the required correlators are obtained as: 

(T T (A n + ik n+1 ) {x, t) (A n - ik n+l ) (0, 0)) 



(6.13) 



[na 



x 2 + (M a T) 2 



x 2 + (u b ry 



Vn,l 



Vn,2 



x cos(Q n x + $ n (r/x)) + 



,2tt 



x + (u a r) 



X 2 + («6t) 2 / 



cos(— + Q' n x + *„(r/x)) 
3a 



(6. 



where n = 1, 4, 6. The exponents are given by: 





1 

2K a 


+ ~ 


^1,2 = 











(6.15) 
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1 

2^ a 


Vl,2 = 


K b 
6 










(6.16) 


^4,1 


= ^6,1 = 


1 

8K a 


+ ^ 
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^4,2 


= ^6,2 = 


3 
8K b 


3^ 6 
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(6.17) 


?74,1 


= = 


1 

8^ a 


+ T 


?74,2 


= %,2 = 


3 

8^6 
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(6.18) 



It can be checkedthat for u a = K a = K b = 1, one recovers the exponents of the isotropic 
SU(3) spin chaint^, namely u U; i + v n ^ = 1 and r\ n ^ + r\ n ,2 = 2/3. One also has: 



(6.19) 



Qi = 0, Q[ = -7r f m - — 

Q4 = = — Qsi Q'4 = ~~ ~y~ = ~Qe- 
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Recall (S z ) = 5/6 + J 2/3mfc. This allows the determination of all incommensurate modes. 
Finally, the functions: 



$1 



$1 



( - 

X 



2 arctan 



T 



X 



UgT 

X 

1 

- arctan 
2 V x 



T 

X 



= 

H — arctan 

2 V x / 



(- 

,x 



U a T \ 

arctan ( — arctan . 

x J \ x J 



fu b r\ 



(6.20) 



All the preceding results are valid only at T = 0. However, it is useful also to calculate 
the correlation functions for T > 0, in particular in order to obtain NMR relaxation rates. 
To obtain the finite temperature Matsubara correlation functions, we can use a conformal 
transformation since we have two decoupled c = 1 conformal field theories. The explicit 
expression of this transformation is: 



/ 2ir{x + iUjT) \ 

x + iUiT — > pui sinn , 

V PUi J 



(6.21) 



where i — a,b. Therefore, to obtainU the finite temperature Matsubara correlation func- 
tions, one has to make the substitution: 



X 2 + (UiT) 2 



arctan 



UiT 



x 



(pUi 



arctan 



cosh 



/ 2ttx 



\f3ui ) 
tan (2vr^) 
tanh [27r-f- 



cos 



2nr 



(6.22) 



With the help of the above results for the spin-spin correlation functions we can evaluate 
the T dependence of the NMR longitudinal relaxation rate T\, 



^ oc lim / /3 rfre^(T T 5 p + (0,r)5;(0,0)) T . 



(6.23) 



We find, 



1 , K. 



oc (a T^ + ^ L ~ 1 + b T^ + ~f +cT 



Ti 



3K a +K b i 3 _j_ 3 3K b + K a 

24 T 8K ( - Too-. x i J rn sir. T o ^. T o 1 



a ~ SK b 1 _)_ d J* 8K a ~ SK b 



The low temperature exponent is the smallest of the the four exponents above. 



4- Comparison with a spin-1 chain with biquadratic coupling 

In the case of a bilinear biquadratic spin-1 chain defined in Eq. (|3.9|), close to the 
Uimin-Lai-Sutherland point (/3 ~ 1), a mapping onto an anisotropic SU(3) spin chain is 
also possible^. However, there are important differences. First, the expression of the spin 
operators in terms of the Gell-mann matrices is different from the ones obtained in the 
spin-tube case. For the spin-1 case, one has: 
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S* = Xl (6.24) 



These expressions should be contrasted with Eqs. (|2.8|) -( |2.9|) . Although the expressions 
( |6.24| ) lead to incommensurate modes, the expression of the correlation functions is different 
from the case of the spin tube. Second, the expression of the Hamiltonian in terms of A 
matrices in the spin-1 case is different from the expression (|2.17|) . Namely, the Hamiltonian 
r9l) can be rewritten in terms of Gell-Mann matrices as: 



H 



E 



tt(A^ +1 + + A 2 A 2 +1 ) + (1 - ~)Af A? +1 



+ Q^^i+i + Ai Af +1 + + AjAj +1 ) H 2~^(Ai A^+i + A? A* +1 - AfAj +1 - AjAf +1 , 



(6.25) 



Finally, for j3 < 1 the spin-1 bilinear biquadratic chain has a gap and the two component 
Luttinger liquid can only be observed for a large enough applied magnetic field. 

Nevertheless, the two problems have in common the presence of a gapless two compo- 
nent Luttinger liquid ground stated, and the formation of incommensurate modes under a 
magnetic field, so that loosely speaking they belong to the same universality class. This can 
be understood as a consequence of the fact that both models can be related to anisotropic 
SU(3) spin chains. One should note that the formation of incommensurate modes in the 
presence of the magnetic field in the spin-tube is not related to the presence of gapped 
incommensurate modes in the bilinear-biquadratic spin-1 chain!!. In the latter case, the 
incommensurate modes originate from the fact that in the absence of the biquadratic chain, 
the (gapped) modes of the spin-1 chain are at q = and q = it whereas at the ULS point 
the (gapless) modes are at q = and q = 2ti/3. The presence of gapped incommensurate 
modes between these two limits is merely a consequence of the continuity of the transition 
between the Haldane gap phase and the gapless phase beyond the ULS point. On the other 
hand, in the presence of the magnetic field, the gapless modes of the spin-tube or those of 
the spin-1 chain simply move away from 2n/3 similarly to what happens in a single spin-1 /2 
chain. 



B. Is there a magnetization plateau for (S z ) = 1 ? 
1. The Umklapp terms and quantization condition on the magnetization 

In the presence of a magnetic field, one of the central issues is the quantization condition 
on the total magnetization (S z ) for the appearance of plateaus. This condition may be 
investigated by looking at the bosonized expression for the spin-operators ( A.12 ). After using 



the transformation ( |A7|) to take into account a non-zero magnetization, one can rederive 
an expression for the non-SU(3) symmetric perturbations. Contrary to the case of zero 
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magnetization, we cannot assume a priori that the "4fcp" terms are highly oscillating since 
the phase e AtkpX may be compensated by a phase arising from the transformation (|A7|) . This 
can be interpreted as a condition for Umklapp processes between the three different fermion 
flavors. A systematic investigation indicates that the only possible Umklapp terms originate 
from the term Af Af +1 or Af Af +1 in the Hamiltonian ( f4.1| ). These Umklapp terms are: 



cos 



cos 




a + 2i/^0b + a 2 ) 



cos(4W^0 b + a 3 ) 



The conditions that must be fulfilled so that these Umklapp terms be present are: 



<A d > 



3v (A 8 ) _ 1 



x/3 3 
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3' 



for the term ( 3.26 ), 



(A-> 



3x , (A 8 ) 1 2 



x/3 3 
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for the term ( |6.271 ), and 



(A 8 ) 1 1 
— =- = — or -. 

6 3 



(6.26) 
(6.27) 
(6.28) 



(6.29) 



(6.30) 



(6.31) 



for the term ( |6.28| ). If we take into account the fact that (A3) = 0, there are only two 
nontrivial conditions to obtain an Umklapp. These are the condition (|6.30|) that reduces 
to (S z ) = 7/6, and the condition ( |6.31| ), reducing to (S z ) = 1. At these values of the 
magnetization the operators (respectively) ( |6.27| ) or (|6.28|) can in principle open a gap in 
the excitations of the system leading to a plateau in the magnetization curve. In agreement 
with the LSM theorem, this implies a ground state that breaks translation invariance. For 
this to happen, these operators must be relevant, namely, 



for the term ( |6.27| ), and 



{K a + 9K b /4) < 1 



K b < 3/4 



(6.32) 



(6.33) 



for the term (|6.28j ). If one of these conditions is satisfied, a gap exists at (S z ) = 5/6 more 
likely than at (S z ) = 1 or (S z ) = 7/6. 
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2. A theory of the possible plateau at (S z ) = 1 



There seems to be evidence for an extra plateau at (S z ) = 1 in the magnetization curve 
of a 36 sites system at J±/ J = 3 obtained by DMRG in Ref. [18] (see Figure 10 of Ref. p~8| ) . 
However, this small plateau observed at (S z ) = 1 in the magnetization curve may also be 
ascribed to the small system sizefi Here, we will investigate in more details the possibility 
of such a plateau in the framework of the bosonized theory. In particular, we will try to give 
a description of the behavior of correlation functions in the system. 

Due to the presence of the Umklapp term (|6.28 ) the bosonized Hamiltonian describing 
the low energy excitations of the spin tube at (S z ) = 1 is: 



H = H a + H b 

dx 



dx 
dx 



2tt 
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20i 



u b K b (^U b ) 2 + ^(d x ^ 2 



+ 



(27ra) 2 
2u'a 



dx cos v 8 




i , /2, 2tt 
'/•'•cos j l\ / -<p b + — 



(6.34) 



If there is indeed a magnetization plateau at (S z ) = 1, a gap opens in the excitation spectrum 
of the field <p b . However, the Hamiltonian contains no terms coupling </> a and <j) b , so that (f) a 
could remain ungapped. This would lead to a single component Luttinger liquid behavior on 
this plateau, and a power law decay of some spin-spin correlation functions. Since v! < 0, 
the Umklapp term would impose \f\{4>b) — — ■ If this is so, careful treatment of the 
expressions (|A12|) of the bosonized forms of the A operators is necessary when we have 
a gap, on the (S z ) = 1 plateau. To eliminate gapped states completely, one can use the 
following expression (see appendix . 



a; 



3Ci,3 c i,3j 



(6.35) 



for S z instead of Q4.12|) here we use the constraint ( f4.11| ). This equation indicates we have 



only the gapful excitations (j) b in A (x). One has 



(A x +zA 2 )(x) 
A s (x) 



Where: 
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(6.36) 



c 2 



(e'V i(<t>b-{<t>b))\ 



(e %: 



i2<J U<h-{<h)) 



(6.37) 



We do not give the expressions for the other operators since they show exponential decay ex- 
cept for A 3 , never enter spin correlation function. We find then that translational symmetry 
is broken on the (S z ) = 1 plateau, with a period of 4 for the ground state. 
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This is in agreement with the LSM theorem that rules out a translationally invariant 
plateau at (S z ) = 1. We also see that an additional period of 4 will appear in the correlation 
functions showing a power law decay. The correlation functions for the operators of Eq. 
is: 



(T r (A 1 + zA 2 )(a;,r)(A 1 -zA 2 )(a; , ,0)) 



(na) 
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Using the expressions (|2.8|) and ( |2.9|) of the spins in terms of A matrices we find that the 
correlations (S+S~,) show an exponential decay whereas the correlations (S* p S*t p ) follow a 
power law decay. We have the following expressions for the equal time spin-spin correlation 
functions: 



C 



7T 



TO 



9 



(6.38) 



Comparing with figures (6) and (7) of Ref. 18], one sees that such behavior is not ob- 
tained in numerical calculations. This leaves two options: one is that the system size (36 



sites) in Ref. [18] is too small to observe the finite but large correlation length. This is not 
unreasonable, since Kb could be only slightly smaller than 3/4. The second possibility is 
that the plateau at (S z ) = 1 is an artifact of the small system size. In Sec. |V1 B 3| it will 
be shown using DMRG for systems of up to 120 sites that it is the latter possibility that is 
obtained. 

If we wish to obtain non-trivial plateaus smaller values of Luttinger parameters K a or K^ 
are needed. This could be achieved by adding a sufficiently strong antiferromagnetic Ising 
term along the chain, 



or an extra coupling, 



N 3 

E\ QZ QZ 

/ < D i,p D i + l, f 

i=l p=l 



N 3 3 

X] X! Si,pSi+l,q- 
i=l p=l q=l 



(6.39) 



(6.40) 



The extra plateaus would lie at (S z ) = 5/6, 1, 7/6 and are allowed by an extended LSM 
theorem in the case of a periodic ground state. 
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3. Zero gap for (S z ) = 1 and finite size scaling of DMRG results 

We continue the density matrix renormalization group3i3 study of the three leg ladder 
at (S z ) = 1 given by Tandon et aZ.EH. Their results for finite chain length show there might 
be a plateau at (S z ) = 1 but the system size is too small to draw definitive conclusions^. In 
this section, we show that the apparent plateau at (S z ) = 1 is indeed a finite size artifact. In 
finite size study, there is a finite energy gap between any two nondegenerate energy levels. We 
therefore use 1/N scaling to show that the energy gap scales to zero in the thermodynamic 
limit for low energy excitations with AS tot,z = and AS tot,z = 1, respectively. 

The system is not dimerized since there is no gap in both AS tot ' z = and AS tot,z = 1 ex- 
citations. Besides the information that there is a certain kind of gapless excitations obtained 
from LSM, numerical analysis gives information on whether AS tot,z = and AS tot,z = 1 
excitations are gapped, and shows whether the system is dimerized. We conclude from both 
the DMRG calculation and the RG analysis that AS tot,z = 1 excitations are gapless for J± 
much smaller than J to J± bigger than J and there is no magnetization plateau at (S z ) = 1. 

We use periodic boundary condition for the rung and open boundary condition for the 
other direction, the i direction in Eq.( |2.1| ), in our DMRG calculation. We calculate the 
lowest energy states for S z = N — 2, N — 1, N, N + 1, N + 2 and denote the lowest energy 
in each S z sector as E (S Z ) and the second lowest energy as Ei(S z ). We have set h = 
in Eq. (|2.1|) and we use J±/ J = 3 to calculate E (S Z ) and Ei(S z ). In DMRG we keep 200 
optimized states and biggest truncation error is 10 -6 . We make calculations for even chain 
length from iV = 4 to iV = 40. 

For each length N we calculate the AS tot ' z = gap by Ei(S z )-E (S z ), and the AS tot ' z = 
1 gap by E (S Z — 1) + E (S Z + 1) — 2E (S Z ). For m — 1, we have plot Ei(S z ) — E (S Z ) vis 
1/N and E (S Z - 1) + E (S Z + 1) - 2E (S Z ) vis 1/N at magnetization S z = N in Fig.|. The 
fittings to the second order of 1/N in the figure show that the gaps scale to zero in form 
Ex{S z ) - E (S Z ) ~ 1/N and E (S Z - 1) + E (S Z + 1) - 2E (S Z ) ~ 1/N when iV -> oo. 

We can see the zero gap excitations also by analyzing the spectrum. Around (S z ) = 1 at 
length iV and S z we will have 3N/2 — S z doublets given in Eq.( [2.3|) . When 3iV/2 — S z is odd 
the ground state is doubly degenerate due to the permutation symmetry of the two kinds of 
the doublets and the two ground states have parity — respect to the i to iV + 1 — i reflection 
symmetry. When 3iV/2 — S z is even, the ground state is unique, but such unique ground 
states for iV and for N + 4 have different reflection parity. It suggests that there is no gap. 
There is no such a parity change from iV to iV + 4 for ground states of gapped translational 
invariant systems or dimerized system. It supports the basic picture in previous sections: 
when we increase or decrease S z , we put in or take out gapless quasiparticles (doublets here), 
these gapless quasiparticles have different parity on different energy levels. These parity and 
degeneracy can be obtained by further detailed analysis of the low energy excitations of the 
Luttinger liquid.0 

We conclude that there is no plateau at (S z ) = 1 since there is no gap for AS tot,z = 
1 excitations. To obtain numerically the non-trivial plateaus, we need smaller value of 
Luttinger parameters K a or K^. If we add a sufficiently strong Ising term along the chain 
direction (in J), the condition is satisfied and we can obtain such plateaus at (S z ) = 5/6, 
1, and 7/6 as we predicted out in previous subsections. These plateaus are allowed by an 
extended LSM theorem in the case of periodic ground state with a certain periodicity, as 
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explained in section III. DMRG calculation for (S z ) =5/6 and 7/6 asks for more technique 
effort!. Calculating numerically the Luttinger liquid exponents^ is also important to fulfill 
the understanding of trileg ladder in future studies. 

VII. CONCLUSIONS 

In this paper, we have analyzed the strong coupling limit of the three chain system with 
periodic boundary conditions in the presence of a magnetic field, using a mapping on an 
anisotropic SU(3) chain. A straightforward extension of the LSM theorem allowed us to 
locate the possible magnetization plateaus. Then, we applied bosonization and renormal- 
ization group techniques to show that for 1/2 < m < 3/2, the system would be described 
by a two component Luttinger liquid. This allowed us to obtain the spin-spin correlation 
functions of the system, and to follow the position of the various incommensurate modes 
in the spin spin-correlation function as a function of the magnetization. Finally, we have 
predicted the temperature dependence of the NMR relaxation rate in this region. We have 
also discussed the relation of the two component Luttinger liquid we obtained with the two 
component Luttinger liquid phase of the bilinear biquadratic spin-1 chain, and shown that 
despite the similarity of the two problems, the two phases are not related to each other. 
We considered the evidence for a magnetization plateau at m — 1 in the framework of our 
bosonized description, and concluded that if such a plateau exists, the ground state at m = 1 
should break translational symmetry with a period of two lattice spacing. We obtained the 
correlation functions of such a ground state as well as the average value of the magnetization 
at each site, and found that if such a plateau does exist the ground state would exhibit some 
kind of antiferromagnetic order. The numerical simulation of Ref. [18] shows no evidence for 
such antiferromagnetic order thus pointing to an absence of plateau at m — 1. To clarify 
whether or not there is a plateau at m — 1, we calculate energy gap around m = 1 using 
DMRG method and fit the data in the linear function of the inverse system size. No gap was 
found in the thermodynamic limit, therefore, we conclude that there is no plateau at m = 1. 
One obvious direction to extend our work is to calculate numerically the Luttinger liquid 
exponents of the three chain ladder with periodic boundary conditions. In appendix [D], 
indications can be found on how these exponents could in principle be extracted. The study 
of these exponents in the case of anisotropic integrable SU(3) spin chains would also be 
interesting. The present paper has been almost entirely concerned with the limit J± ^> J, 
and it would be interesting to investigate whether the behavior we have obtained is valid 
also in the opposite limit, J± <C J. It would also be interesting to study anisotropic gen- 
eralizations of the spin-tube model and check for plateaus at m = 1 as well as extend the 
analysis to higher n spin-tubes. Finally, the relation between the model we have considered 
and classical statistical systems such as Coulomb gases may be worth analyzing in more 
details. 
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APPENDIX A: ABELIAN BOSONIZATION OF THE SU(3) HUBBARD MODEL 



In this section, we apply abelian bosonization to the SU(3) Hubbard model: 

H h = -*X][cJ,n C i+l,n + h.C.] + — n i,nni,m- (Al) 

i,n i,n^m 

In the strong-coupling limit, the Hubbard Hamiltonian with one fermion per site projected 
onto the low-energy states becomes simply the Heisenberg Hamiltonian as can be seenEl 
by considering perturbation theory in the hopping term for U ^> t. Only second-order 
perturbation theory survives and the effective Heisenberg coupling is J = t 2 /U. 

In the continuum limit, in terms of the right and left moving fermions introduced in Sec. 
1V| , the free Hamiltonian H t can be rewritten as: 

H t = -W / dxJ2(' l pR,ndxi ! R,n ~ ^L,n d x^Ln), ( A2 ) 

where v = 2ta sin(A;^a) is the Fermi velocity. In the following, we are working at a filling of 
one fermions per site. This implies k F = ?r/(3a) and vf = \/3ta. 

Using the standard dictionary of abelian bosonizationEa we express ipR(L) n m terms of 
the bose fields 4> n and their duals 9 n , for each flavor n = 1, 2, 3, 

^ Rn (x) = - * eW'O+M*))^ 
y/Zira 

^ Ln (x) = -^e^W-^W)^, (A3) 

where ?7_R(L) n are the Klein factors ensuring the proper anticommutation relations among 
fermion operators^. One has: 7rn n (x) = d x 6 n and [(j) n (x), H m (x')] = i5 n ^ m 5(x — x'). 
The non interacting Hamiltonian is straightforwardly rewritten as: 

Ht = Y.v fp [(vrn n ) 2 + [dAnf] (A4) 

n=l J 

And the fermion densities as: 

p n (x) = + - e 2 ^ + 6 e- 2 ^ (A5) 

7r 2ira 2ira 

where k F = 

The Hubbard interaction, V = (Ua/2) J2n^m I dxp n (x)p m (x), is rewritten in terms of the 
fields 4> n , 

f U 2U 

V = dxV] -^d x (j) m d x (j) m + — — — cos 2(<j) n - <f> m ) (A6) 

Instead of working with fields 0i, fa, 03 it is convenient!^! to introduce the transformation: 
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and similarly for the conjugate fields. The field (j) c describes the charge excitations, whereas 
the fields cj) a ,b describe the SU(3) spin excitations. We recover in particular the fact that the 
SU(3) spin excitations are described by a conformal field theory with C = 2 whereas the 
charge excitations have C = 1. The charge and spin sectors of the Hubbard Hamiltonian 
are then completely separated, 



H = H c + H, 
dx 
~2n~ 



H, 

H, 



where 



u, 



i=a,b 

2g 
+ (2na) 2 



(A8) 
(A9) 



dx 



cos(\/80 a ) + cos V2((f) a + V3(f) b ) + cos v2(0 a - \/306 



(A10) 



U r Kr. 



U.K. 
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v F + 



v F - 



2Ua 

71 

Ua 

7T 



g = Ua 



(All) 



Note that the Hamiltonian describing the charge modes contains no Umklapp term that 
could lead to a gap opening. This is due to the fact that Eq. (|A5|) has been truncated at the 
2k F harmonics. In a more complete expression!^, higher harmonics would appear and would 
give a higher order 6fcp-Umklapp. This Umklapp term can also be derived by perturbation 
theory@. This Umklapp term is of the form cos2\/30 c and is irrelevant for U/t 1. This is 
confirmed by numerical simulations^! which show that the charge gap in a SU(3) Hubbard 
model open only for U > 2.2t. When considering RG equations (see Sec. [IV]) we shall see 
that in the spin sector, for U initially positive, K s renormalizes to 1 and g renormalizes to 
0. This implies that the spin sector of the SU(3) Hubbard is described by a c = 2 conformal 
field theory perturbed by a marginally irrelevant operatoilllii. 

Of course, we also need a bosonized expression for the SU(3) spin operators. This 
can be derived from the continuum limit of the definition ( 4.12j) of these operators, recall: 
A a (x) ~ a _1 A", x = ia (a = 1 • • • 8). We obtain: 



A 1 
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cos V29 a 


2 cos V2d> a + e^ x e~ 2t ^ + e~ l ^ x 


2i^' 
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sin \pld a 


1— 2tt 9, 2ir 9, ^b ' 
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A 4 (x 
A 5 {x 
A 6 (x 
A 7 (x 
A 8 (x 




2 cos 



2 cos 



2 cos 



2 cos 




( ±<L 
\V2 



„ ,( $a.-*k\ 



4>a I 4>b 



+ e z 3 



+ h.c. 
+ h.c. 
+ h.c. 
+ h.c. 



cos 



+ h . c. 



(A12) 



Where A a (x) = for x = na. Using these expressions, one can derive immediately the 
expressions ( |4.20|) -( [4.21 ). In the limit U — > oo, one must note that the expression of Ag(x) 
has to be modified. The reason is the following: for U — > oo, one has C1C1 + C2C2 + C3C3 = 1 on 
each site. As a result, A 8 = (cjci + c^ — 2c3C 3 )/v / 3 can be rewritten as: A 8 = (1 — 3c\c 3 )/\i r 3. 
Using bosonized expressions, one obtains: 



A 8 (x) = -—d x <p b - 



4>b 



+ H. c. 



(A13) 



Thus, the terms containing cos v^20a drops from the expression of A 8 (x) in the limit U — > oo. 
This means that the SU(3) Hubbard model with a finite charge gap and the SU(3) spin chain 
should have in general different correlations for A 8 . It should be noted that this difference 
should not appear at the isotropic point, where the exponents of the correlation functions 
are identical. However, it is obtained for models in which the SU(3) rotation symmetry is 
broken. 



APPENDIX B: OPERATOR PRODUCT EXPANSION OF MARGINAL 

OPERATORS 

In this section, our aim is to derive the OPE for operators of the form cos(v / 8«.0), and 
deduce the renormalization group equations. 

Let us first recall briefly how operator product expansions can be used to obtain one 
loop renormalization group expansions^. Assume we are given a set of operators closed 
under the OPE 

$ i {x,T)$ j (0)~Y,c ij k (x,T)$ k (0), (Bl) 

k 



in the sense that expression (B 1|) , when inserted in any correlation function, gives the correct 
leading asymptotics for (x, r) — > 0. Denoting [<J>j] the scaling dimension of the operator <3>j, 
defined by : 
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($,(x,0)$ 4 (0,0)> cc I-) 2 ™ (B2) 



x 



then c ijk (x, 0) ~ const. x^H^H*]]. 
Perturbing the Hamiltonian by 

Hi = J2 dxdTg k $ k (x,r) (B3) 
k J 

One can deduce the one loop renormalization group beta functions directly from the OPE 
( |B1~1) . These one loop renormalization group equations are: 

^ = p k (g) = (2 - [S> k ])g h -7r£Cy*ftfc (B4) 



Where we define 



C ii = a 2 [ ^-c ijk (a cos 9, - sin 0) (B5) 

JO 27T U 



The Equations [BJ are a slight generalizations of those that can be found in Ref. |65|, in which 
we have allowed for a function Cij k (x,r) that depends both on x 2 + u 2 t 2 and urjx. The 
set of operators has to be closed under the operator product expansion (i.e. they have 
to form a closed algebra), otherwise new operators would be generated under RG, and new 
OPEs would have to be derived. We thus need to retain the smallest closed algebra that 
contains all the operators that appear in our problem. In our case, we have to retain the 
operators cos\/8c?i-0 as well as the operators (d x 4> a ,b) 2 - 

In order to derive the OPE for the cos\/8aj.</> operators, we use the following identity: 

j<x4 = : e w : (B6) 

Where : . . . : represents normal ordering. This identity implies: 

e is/83,$(x,r) e -iV8a.$(0,0) _ e -4((a.(^,r)-^(0,0))) 2 > . e tV8a.($(x,r)~$(0,0)) . 

We have: 

e -4((a.(0(x,r)-0(O,O))) 2 > = ( ^ \ R /ggN 

\x 2 + (ut) 2 J 

And we can expand the normal ordered product, yielding: 

. el Vsd.($(x,r)-$(o,o)) ._ l _ 4(^(^0(0, o) + r<9 T 0(O, 0))) 2 (B9) 

This leads to the OPE Q. 
Now consider: 

e lV8S-$(x,r) e iVs0.$(O,O) _ e -i{(d.$(x,r)+0.$(O,O))) 2 ) . e iVS(a.$(x,r)+0.$(O,O)) . (BIO) 



32 



Rewriting it 

((a.fe r) + /3.0(O, 0))) 2 ) = (((a + /3).0) 2 > - 25/3((0 2 (O, 0) - 0(x, r)0(O, 0))), (Bll) 



we have 



/ \ -2KS.0 

c iV8a.$(x,r) c iV80.$(O,O) = | Q j e %y/%{5+0)4 (B12) 



^Ac 2 + w 2 r : 

The OPE for the cosines can be obtained trivially. There is no need to calculate the OPEs 
of the operators {d x cf) a ^) 2 with the operators cos>/8<3j.0. It can be shown that these OPEs 
only reflect the dependence of the scaling dimensions of the cosines with K a ^. Therefore, 
we have obtained all the OPEs needed to derive the renormalization group equations. It 



is then a simple matter to write the renormalization group equation using the formula [B4 
Following the procedure described in Ref.0, one obtains: 

— = i2 -2K a ) yi - — 

dy 2 ( 2 _Ka_ z K,\ vm 



dl V 2 2 J a 2 
dl \ 2 2 J yi 2 
dl \K a ) 2 8 8 

d f 1 \ -_ 3 yl + 3 yl (B13) 



dl \K, 



where yi = A few remarks on these equations have to be made. In Ref. |65|, the OPE 
depends only on the distance between points. In our case, the OPE also depend on the angle 
between the segment joining the points and the horizontal axis. Since in the derivation of the 
RG equations one integrates over the ring a < r < ae dl the angular part of the integration 
cancels the terms d x <pd T <p and gives a n/2 factor for the terms (d XjT (j)) 2 . The second important 
remark is that in our equations, we are working with 2/2(0) = 2/3(0). It can be checked that 
this condition is preserved by the RG flow and that under such condition no terms d x <j) a d x <f)b 
are generated. Finally, if we expand for small 2/4,2/5 , we have K a = 1—2/4 and Kb = 1 — 2/5. 
Putting this in Eqs. (|B13| ) we get the renormalization group equations 



APPENDIX C: RENORMALIZATION GROUP EQUATIONS IN PRESENCE OF 

THE EXTERNAL MAGNETIC FIELD 

In the present section, we want to extend the derivation of the renormalization group 
equations of appendix |B| to the case of a non zero effective magnetic field. As explained in 
section [VI|, it is convenient to first perform a Legendre transformation and work at fixed 
magnetization. Then, the RG equation for the magnetization becomes trivial, but the 
RG equation for the magnetic field is not. Similarly to the zero magnetic field case, the 
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renormalization group equations can be obtained via operator product expansion. The only 
difficulty is that the problem is not a priori translationally invariant. 

The relevant operator product expansions are obtained by the method of appendix [B|. 
One has: 



cos 



x, r) + ux)) cos (\/8a.((j)(x', r) + ux'] 



{x 2 + U 2 T 2 ) 2 

x (— v^sin (V8(u.a)x) (xd x (a.(f) — rd T (a.(f)) — 2 (a.(xd x (j) + rd T (p)) cos (Vs(u.a)x) j (CI) 



In our case, one must take u = — 7rm b (0, 1). The term in rd T (a.(j)) disappears upon 
angular integration. On the other hand, the term in xd x (a.<j)) leads to a renormalization of 
the applied magnetic field. The angular integrations lead in general to Bessel Functions. 

The second useful OPE is: 



-2Ka.f3 



e iV8a.($(x,r)+ux) e iV80-(<p(x' ,0)+ux') 



X — X') 2 + U 2 (t — T /N|2 



,tV8(a+0) . {$( ,0)+52±2- ) 



(C2) 



These OPEs allow us to deduce the renormalization group equations for K a , K b , yi, y 2 , h 
in the form: 



d / 1 
Jl \K~ a 
d / 1 
Jl \K~ b 



1 2 1 2 t / /r\ V^3 . 

2^1 + 4?/2^o(vrm b (/)— a) 

3 fs 
-yl J (7rm b J -a(l)) 



% = (2 - 2K a ) Vl - ly 2 2 J (7im b (l)V3a) 



dl 

dm 
dl 

dh 
111 



1 3 \ 
2 - -K a - -K b J y 2 - ymJ^—m b (l)a) 



— = J —ylJi(7i\/6m b (l)a) 



(C3) 



APPENDIX D: DETERMINATION OF THE EXPONENTS OF THE 

BOSONIZED HAMILTONIAN 

In this section, we discuss the determination of the exponents of the spin-tube. We 
have shown previously that in weak coupling the model flows to a two-component Luttinger 
liquid fixed point. In strong coupling some alternative techniques are needed to determine 
the Luttinger liquid exponents from thermodynamic quantities. Note that in the isotropic 
SU(N) Hubbard model case^l with a charge gap, one needs only the spin velocity since 
the spin exponents are constrained by SU(N) invariance. Here, SU(3) symmetry is broken 



34 



and we expect two different velocitiesw a , Ub and two exponents K a ,Kb so that we need four 
independant quantities. 

Suppose that we have a two component Luttinger liquid described by the Hamiltonian: 



H= V [*L 

i=a,b 



(Dl) 



In the case of the spin tube, we can rule out terms of the form n a Ilft and d x <p a d x (j)b since 
we know that they are not present in the bare Hamiltonian and not generated upon RG. 
Moreover, we know that X38 — which guarantees the absence of terms d x <p a d x (f)b in the 
Hamiltonian. 

The usual technique to determine the Luttinger liquid exponents is to consider the energy 
change induced by taking 7i(Yl a ^) = (p a ,b/L. In terms of the Luttinger liquid parameters, 
this energy change is given by: 

This energy change is related to the change of ground state energy caused by taking twisted 
boundary conditions. Let us discuss in more details these twisted boundary conditions in 
the specific case of the spin tube. Using the bosonization formulas, one sees easily that: 

A 1 + ik 2 oc e lV29a 
A 4 + *A 5 oce 

A 6 + zA 7 oce l (~^' V3 "°J (133) 



Therefore, imposing (7rn a ) = ip a / L and (nllb) = <fb/L amounts to imposing the boundary 
conditions: 

(A 1 + ik 2 )(L) = (A 1 + ik 2 )(0)e lV2 ^ 

(A 4 + iA 5 )(L) = (A 4 + iA 5 )(0) e ^ a/V2+ ^ b) 

(A 6 + ik 1 ) (L) = (A 6 + ik 7 ) (oy'Vlw-WVS) (D4) 

As an aside, one should remark that the transformation:n(x) — > U(x) — f(x) is realized 
by the operator U = exp (— % J dxf(x)(p(x)). This operator can also be written as: U = 
exp (i J dxF(x)d x (p(x)) where / = Twisted boundary conditions correspond to f(x) = 



a/L. Therefore, an operator generating states satisfying boundary conditions (P4j) acting 
on states satisfying periodic boundary conditions can be built in the continuum. A lattice 
version is easily constructed, giving an operator of the form: 



L (n-1 



U(ip) = exp -, £ ^(yv^ + yvV*) (D5) 



n=l 



One can check that these lattice operators acting on states that satisfy periodic boundary 
conditions generate states that satisfy the boundary conditions ( |D4| ) directly on the lattice. 
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This guarantees the existence of states satisfying the twisted boundary conditions ([DJ). The 
generalization of this construction to the SU(N) case is trivial. Instead of A 3 ' 8 one has to 
consider the Maximal Abelian Subalgebra (MASA) and builds the operators corresponding 
to U(ip). In the case of the spin-tube, the energy change of the ground state obeys the 
condition: 

LSE(ip a , cp b ) = L5E(<p a , 0) + L5E(0, <p b ) + o((p 2 a , cp 2 b ) (D6) 

In order to determine completely u a , Ub, K a , Kb suppose that one places the anisotropic 
SU(3) chain in fields that couple to the A 3 and A 8 components of the spin, 

H = E E ^ A? +1 - h 3 E A, 3 - h 8 E A, 8 (D7) 

i a i i 

Then, the Hamiltonian becomes in the continuum: 

H = H + E (D8) 



7T 



where h a = h 3 , hb = h 8 . Then, one has: 

71 U v 

Thus, one has: 



K (D9) 



(A 3 ) = K a /u a h 3 

(A 8 ) = Kb/u b h s (D10) 

This is sufficient to extract the parameters of the two component Luttinger liquid asso- 
ciated with the anisotropic SU(3) spin chain. However, we started from three coupled spin 
1/2 chains with periodic boundary conditions. To extract the two-component Luttinger 
liquid exponents for this problem, we need to express the preceding formulas in terms of the 
original spins. Re-expressing the twisted boundary conditions in terms of the original spins 
is elementary if one remembers that when we choose: S z = (5/6 — A 8 /v / 3), we have: 

(a 6 - i\ 7 )i = J2f~ ls £ P 

(A 4 -^)„ = A^ j2 ( P -D 5 ,+ 

{\ 1 -i\% = - l -Y,f (p - 1) si P (dii) 



The following expression for A 3 can also be obtained: 

((S 2 S+ - S 2 S+) + (S 3 S+ - S+S 2 ) + (S^S+ - S+S3)) P s *=i/ 2 , (D12) 



A, ' 
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where Ps z =i/2 is the projector on the subspace S z = 1/2. Physically, A3 is proportional to 
the spin current in the transverse direction. It is therefore +1 for positive chirality and — 1 
for negative chirality. It can also be rewritten as: 



A, = 



(s 2 x s.y + (s 3 x s 2 y + (£1 x s 3 y 



=1/2 



(D13) 



With these expressions, it is in principle possible to obtain numerically the Luttinger liq- 
uid exponents for a general three leg spin ladder with periodic boundary conditions under 
magnetic field between the (S z ) = 1/2 and (S z ) = 3/2 plateaus. 
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TABLES 

TABLE I. Correspondance between the notations of the present paper and those of K. Tandon 
et al. 





States 








Present paper 


|1> 


|2> 


13) 




K. Tandon et al. 


3 1 7') 


f 1 5') 


|1> 






Operators 








Present paper 


n 


n 




T z + 1/3 


K. Tandon et al. 




JL- 


fR- 
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FIGURES 
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FIG. 2. The energy levels of a single triangle as a function of the magnetic field. Solid lines 
correspond to states with S=3/2, dashed lines to states with S=l/2. One observes the level crossing 
between the state with S z = 3/2 and the states with S = 1/2, S z = 1/2 as the magnetic field is 
increased. 
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FIG. 3. The renormalization group flow with initial conditions (5.3) and yo = 0. The couplings 
constants <7i,ff2j<?3 are renormalized to 0, whereas 54 — > g\, and 55 — > g?-. The system therefore 
flows to a two component Luttinger liquid fixed point. 
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FIG. 4. The renormalization group flow with initial conditions and j/q = 1. The couplings 
constants gi,g2,93 are renormalized to 0, whereas 54 — > <?|, and 55 — > 55. The presence of a 
marginal perturbation preserving SU(3) symmetry does not suppress the two component Luttinger 
liquid behavior. However, by comparing with Fig. ^, it is seen that it changes the exponents at 
the fixed point. 
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Fig.5 Citro et al, PRB 
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FIG. 5. Gap scalings for magnetization m = 1. Lowest two energies Eq(S z ) and Ei(S, 
each S z obtained by DMRG are plotted into AS tot ' z = gap scaling: lOp^S^) - E^Sz 
1/L), and AS tot > z = 1 gap scaling: E (S Z - 1) + E (S Z + 1) - 2L (S Z ) vis l/L with S 2 
We have magnified the AS tot ' z = gap by times to make the figure clear. The linear fittings are 
- E (S Z ) ~ 1/L and - 1) + £b(S z + 1) - 2L (S' Z ) ~ 1/L in thermodynamic limit. 

The chain lengths calculated by DMRG are L = 8, 12, . . ., 40. 
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